The Swift-Hohenberg equation requires 
non-local modifications to model spatial 
■ pattern evolution of physical problems 
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Abstract 



> 

I argue that "good" mathematical models of spatio-temporal dynam- 
ics in two-dimensions require non-local operators in the nonlinear terms. 
Consequently, the often used Swift-Hohenberg equation requires modi- 
fication as it is purely local. My aim here is to provoke more critical 
examination of the rationale for using the Swift-Hohenberg equations as a 
' reliable model of the spatial pattern evolution in specific physical systems. 

in • 1 Introduction 



Consider the spatio-temporal dynamics of systems with a very large horizontal 
O extent, when compared to their height. For two definite examples, I will refer 

to Rayleigh-Benard convection, 



X| ■ ^(f +u-Vu) = -Vp> + Ra8e z + V 2 u, 

|f + u - ve = w + V 2 6>, U 

where Pr and Ra are the Prandtl and Rayleigh numbers respectively, and also 
refer to a toy set of partial differential equations, 



<)a 



(1 + V 2 ) a-ab, 



§ = - (a - a^ 2 + a 2 V 4 ) b + a 2 . 



For values of the parameters near some critical value, a continuum of modes, 
an annulus in Fourier space |k| w const, become linearly unstable. Through 
physical nonlinearities these modes not only saturate, but also interact. The 
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interaction of the continuum of critical modes is the main feature of interest as 
it determines the spatial patterns in the horizontal. In order to explore pattern 
evolution, many researchers IE, |7[ [TOj, [l3| for example] have invoked the 
Swift-Hohenberg |l(| equation, 

d -A = , L A-(kl + ^) 2 A- lA \ (3) 

as a useful system to study. Whereas it is indeed instructive to examine the 
Swift-Hohenberg equation, I contend that the Swift-Hohenberg equation is de- 
ficient as an accurate and reliable model of most specific physical planform 
evolution problems. 

Instead, I propose that it is generally necessary to incorporate non-local 
nonlinearities into a low-dimensional model of planform evolution. Specifically, 
I recommend an equation of the form 

^=^A-(kl + V 2 ) 2 A-AG*A\ (4) 

where Q* is some radially symmetric convolution. My argument to support 
this recommendation has three facets. Firstly, symmetry considerations permit 
non-local operators in a model (§2). Secondly, non-local operators are naturally 
generated in systematic methods of modelling nonlinear pattern evolution |Q 
(§3). Indeed, both Swift & Hohenberg [|l6| and Bestehorn et al jl| naturally 
encountered non-local terms in their derivations; but they heuristically argued 
to replace them by local nonlinearities. Thirdly, the range of Fourier harmon- 
ics generated by the nonlinearities is fundamentally different in two-dimensions 
than in one-dimension (§4). In the two-dimensions of a planform evolution 
problem nonlinearities generate a continuous disc in Fourier space of harmon- 
ics, whereasin one-dimension the generated harmonics congregate in discrete 
lumps. This difference requires a more sophisticated treatment of the two- 
dimensional planform evolution problem, one that necessarily leads to non-local 
nonlinearities. 



2 Symmetry arguments permit non-local effects 

The rationale behind the use of the Swift-Hohenberg equation (|^) is firstly that 
for parameter values near critical, /j«0, the spectrum of the linear terms, A = 
fjb — ffeg — |k| 2 ) , matches, near the critical wavenumber fcoj the spectrum of the 
physical problem under study. Secondly, that the cubic nonlinearity is typical 
for the symmetry of many systems and is needed to stabilise finite amplitude 
dynamics. The symmetry invoked is that of sign: A — ► —A leaves the Swift- 
Hohenberg equation (|^) unchanged; as does a — > —a in the toy problem (0), and 
(8, w, z) — > (—6, —to, —z) in Rayleigh-Benard convection ([j]). This symmetry is 
only consistent with odd functions. However, this symmetry applies to the field 



2 



as a whole, not to every part of it individually. Thus a wide variety of cubic 
functions are potentially permissible: for example, non-local functions such as 
A(x) 2 A(x + e x ) or 

J K(x, x', x", x"')A(x')A(x")A(x"') dx'dx"dx'" , 

where x = (x, y) etc. 

However, translational and rotational symmetry in space requires that any 
non-local effects must be expressible as radially symmetric convolutions. The 
local nature of typical models of physical dynamics implies that harmonics are 
forced by locally expressible functions. But the feedback from the harmonic to 
the critical modes, essential for stabilisation of finite amplitude dynamics, need 
not be local. Indeed a number of systematic studies have shown that mem- 
ory, either temporal or upstream, are needed in the low-dimensional modelling 
of forced dynamics || or shear dispersion in varying channels [|l5| In a 

planform evolution problem such a "memory" of horizontal structure, occurring 
through fluid convection for example, would manifest itself as the non-local con- 
volution in a cubic nonlinearity, such as the AQ*A 2 term in (^) where Q is some 
radially symmetric kernel. 

3 Non-local nonlinearities are natural 

The evolution of the Swift- Hohcnberg equation contains a thin annulus of critical 
modes, |k| ss fco, near onset. It also contains a wide variety of non-critical, 
exponentially damped modes, |k| 56 fco- Thus the interesting long-term evolution 
of the critical modes is embedded within the equation along with the evolution 
of many uninteresting modes. It is for just such a scenario that I developed 
the concept of an embedded centre manifold [ fl4| . There I show that adiabatic 
iteration, namely the repeated application of adiabatic elimination jll], |l8|, [lTj , 
effectively embeds the critical modes of a slow manifold into the dynamics of a 
higher-dimensional system. In a pattern evolution problem, the state space of 
the higher-dimensional system consists of all the modes in the two-dimensional 
plan, whereas the slow manifold is composed of just the annular neighbourhood 
of the critical modes. 

For the toy problem (Q) with a?, = 1 — a\ = 2 and 012 — 1, adiabatic 
iteration leads to ([|) as a first non-trivial approximation where 

[(l-V 2 ) 2 -r]g = <5(x). (5) 

For example, if r = then Q = j^Ko(x) ★ Kq(x) in two-dimensions and Q = 
i(l + |x|) exp(— \x\) in one-dimension. 

Now the left-hand side of (||) comes directly from the spectrum of the expo- 
nentially decaying branches of the linearised problem. The only way to avoid a 
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Figure 1: schematic diagram in wave- number space of the leading order non- 
linear interactions between modes, ki and k2, on the critical circle |k| = fco 
through the nonlinear generation of the harmonic at point B. 



Figure 2: schematic typical amplitude spectrum of two-dimensional pattern 
evolution where the critical modes with |k| = fco have typical amplitude e. 



non-local operator is if the solution of the analogue of (g), in the given phys- 
ical system, is a delta function, Q oc <5(x). Given the typically elliptic nature 
of dissipative and spatially symmetric operators, this can only happen if the 
operator on the left-hand side of the analogue to (||) is constant as a function of 
wave-number |k| (or more generally, constant on each branch). Generally this 
will not occur; for example, in Rayleigh-Benard convection (Q) with stress-free 
boundaries and Pr = 1 the spectrum of the mth branch is 

X = -™ 2 -\\c\2 + ,f~R^ — J k \ (6) 

^Ikl 2 +m 2 

Thus it is generic that systematic modelling naturally leads to non-local 
nonlinearities such as shown in (^). 



4 Forced harmonics in Fourier space 

Here I argue that the richness of non-local nonlinearities are generally necessary 
for an accurate model of a physical problem. The argument rests on the basic 
nature of the power spectrum in the two-dimensional problem. As shown in 
Figure ^, two critical modes with wavenumbers ki and k2 on the critical circle 
|k| = fco interact through quadratic nonlinearities to generate a harmonic at B, 
wavenumber ki + k2 . In the toy problem (Q) the relevant quadratic nonlinearity 
is the a 2 term in the b equation; in convection it is the advection terms, u • Vu 
and u • V0, on the left-hand sides of (Q). Then such a forced harmonic interacts 
with the — ki critical mode, via the ab term in (|J) or the advection terms in ([!]) 
for example, to generate a forcing of the k2 component in the critical modes. By 
varying ki and k2 independently around the critical circle, every wavenumber 
within the disc |k| < 2fco is independently forcedjj] 

Thus the amplitude spectrum of typical planform evolution evolution is as 
shown in Figure ^[ Observe that this amplitude spectrum is qualitatively dif- 
ferent from that in one-dimension. In one-dimension the amplitude spectrum 
consists of discrete lumps at integer multiples, nk 0l of the critical wavenum- 
ber. For example, this property of the one-dimensional problem is crucial in 

1 lt is symmetry, for example a — > —a, which ensures that two critical modes together do 
not directly force another critical mode. 
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proofs that a Ginzburg-Landau equation is relevant to one-dimensional pattern 
evolution |6|. However, in two-dimensional problems the amplitude spectrum is 
considerably richer. Consequently, in order to model the leading order, physi- 
cal interactions among all these modes, it is necessary to determine the forced 
harmonics over the entire disk |k| < 2fco. Because these harmonics naturally 
decay with a rate which depends upon wavenumber, then in physical space this 
determination can only be done via a non-local convolution. For example, in the 
toy problem we solve (||) for use in (||) in order to account for the variations of 
the decay rate —A = — r + (1 + |k| 2 ) 2 of the harmonics b over the disc |k| < 2ko- 

The best that the Swift-Hohenberg equation (J3[) can do is to approximate 
the functional dependence A(|k|), as above or in (g)^ by a constant. It is easy to 
imagine problems where this approximation would be inadequate. For example, 
if in the toy problem (||) we choose ai = — 4a?2 > 0, then the decay rate of 
the harmonic b has a minimum at wavenumbers |k| = \/2ko and so I expect 
that a square planform would be preferred because the harmonics involved in 
the necessary interactions are not damped as strongly as other stable modes. 
(In convection there can be a minimum in the decay-rate, but for Pr = 1, (^|), 
it occurs for |k| < fco.) The Swift-Hohenberg equation misses such subtleties, 
and as such it cannot be expected to be a reliable model of any given physical 
problem. Instead, non-local terms based on the wavenumber dependence of the 
decay of harmonics need to be used in order to make reliable physical predictions 
from such a model. 

Lastly, I give a further extremely formal argument for (|4|) , one based on the 
new notion of matching centre manifolds in order to systematically develop 
models of specific problems. The idea is to develop a model, of lower-dimension 
than the original dynamical system, whose slow manifold evolution matches 
that of the original system to some order of analysis. This principle is very like 
that of using Pade approximants to sum Taylor series. Here I would propose 
the modified Swift-Hohenberg equation (^) as a model, and then determine the 
necessary constraints on the non-local convolution, in order for the slow 
evolution of the critical modes of (^) to be the same as those of the physical 
problem under study. 

Here the most restricted version of critical modes are those precisely on the 
circle |k| = k$. Thus we express the solution field as, for example, 



A(x, t)= Z^t) exp(ikoO) • x 
Jo 



)# + •••, 

for some complex amplitudes Z. Concentrating only upon the nonlinear inter- 
actions, the evolution of the slow modes would take the form 

3Z /" 27r 

~gf = !~ / (3{4> 1 i>)Z (j ,Z^Z 7!+ ^d'4> 



Recently, Edwards & FauveQ have in essence used a discretised version of of 
this, their equation (8), in studying pattern selection in the Faraday experiment. 
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Indeed their interaction diagram, Figure 4, is essentially the same as Figure [y 
above. In order for the matching to take place to cubic order, it is necessary 
that the interaction coefficient, /?(</>, ip), for the physical problem and for the 
model to be identical. As explained before, this interaction takes place through 
all the forced harmonics in the disc |k| < 2k . In essence, (0, ip) appearing 
in the interaction integral above is a straightforward parameterisation of this 
disc. (For example, rotational symmetry in the problems imply that [3 is purely 
a function of the difference <fi — ip.) Thus, the only way for the model (|J) to 
accurately match a specific physical problem, is for the interaction over all the 
harmonics in the disc to be accurately represented. Thus non-local nonlinearities 
are essential. 

The only freedom this matching principle permits, up to third order in the 
analysis, is the freedom to vary the Fourier transform of Q for wavenumbers 
greater than 2fco. This cannot change the non-local nature of the convolution. 
In general, the Swift-Hohenberg equation requires non-local modification. 
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